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ABSTRACT 



Petri 



fl2005a 



a) 



In this paper we extend the idea suggested previously by 
(paper I and II) that the high frequency quasi-periodic oscillations (HF-QPOs) 
observed in low-mass X-ray binaries (LMXBs) may be explained as a resonant 
oscillation of the accretion disk with a rotating asymmetric background (gravi- 
tational or magnetic) field imposed by the compact object. Here, we apply this 
general idea to black hole binaries. It is assumed that a test particle experiences a 
similar parametric resonance mechanism such as the one described in paper I and 
II but now the resonance is induced by the interaction between a spiral density 
wave in the accretion disk, excited close to the innermost stable circular orbit, 
and vertical epicyclic oscillations. We use the Kerr spacetime geometry to de- 
duce the characteristic frequencies of this test particle. The response of the test 
particle is maximal when the frequency ratio of the two strongest resonances is 
equal to 3:2 as observed in black hole candidates. Finally, applying our model to 
the microquasar GRS 1915+105, we reproduce the correct value of several HF- 
QPOs. Indeed the presence of the 168/113/56/42/28 Hz features in the power 
spectrum time analysis is predicted. Moreover, based only on the two HF-QPO 
frequencies, our model is able to constrain the mass Mbh and angular momen- 
tum obh of the accreting black hole. We show the relation between Mbh and 
ttBH for several black hole binaries. For instance, assuming a black hole weakly 
or mildly rotating, i.e. obh < 0.5G'Mbh/c^, we find that for GRS 1915+105 its 
mass satisfies 13 Mq < Mbh < 20 M©. The same model applied to two other 
well-known BHCs gives for GRO J1655-40 a mass 5M0 < Mbh < TMq and for 
XTE J1550-564 a mass 8M0 < Mbh < 11 Mq. This is consistent with other 
independent estimations of the black hole mass. Finally for H1743-322, we found 
the following bounds, 9M0 < Mbh < 13 M©. 
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Introduction 



High frequency quasi-periodic oscillations (HF-QPOs) are common features to all 
accreting compact objects, be they neutron stars, white dwarfs or black holes. A number 



of recent observations have r e vealed the existence of these HF - QPOs in several 



Dlack hole 



binaries ( Strohmayer 



2001a : 



MacClintock fc Remillard 



2003 



Remillard et al. 



20061 ) ■ 



Usually, a pair of HF-QPOs appears in a 3:2 ratio. If these oscillations are connected to the 
orbital motion of the accretion disk at its inner edge as predicted by several models, these 
QPOs become a useful test of gravity in the strong field regime. 



The 3:2 ratio was first noticed by Abramowicz fc Kluzniak 



fl200lh . In order to 



explain this ratio, they introduced a resonance mechanism between orbital and epicyclic 



motion around 



i eiT bla ck holes therefore leading to an estimate of their mass and spin. 



Kluzniak et al. 



(j2004al ) showed that the twin kHz-QP Os is expla i ned b y a non linear 



resonance in the epicyclic motion of the accretion disk. 



Rebusco 



(12004h developed the 



analytical treatment of these oscillations. This parametric e picyclic resonance model o f 



hydro dynamical modes in the accretion disk was applied by 



Kluzniak fc Abramowicz 



( 120021 ) to some microquasars. They pointed out the 3:2 ratio in the HF-QPOs observed 



in GRO J1655-40, X ' 



binary. 



Strohmayer 



E J1550-564 and GRS1915+105. Moreover, in this latter black hole 



(j2001bl ) reported anoth er pair of frequencies, naraely 69 .2 Hz and 



41.5 Hz, which are in a 5:3 ratio as noticed by|K 
the parametric resonance model. Furthermore 



uzniak &: Abramowicz 



RezzoUa et al. 



(I2OO2I ) and supports 



(I2OO3I ) suggested that the 



HF-QPOs in black hole binaries are related to p-mode oscillation in a non Keplerian torus. 
Nevertheless, the propagation of the emitted photons in curved spacetime can 



also produce some intrinsic 



( ISchnittman fc Bertschinger 



peaks in the Fourier spectrum of the light curves 



2OO4J ). In combination with a verticall y oscillating 



torus, the gravitational lensing effect can also reproduce the 3:2 ratio ( .Bursa et al. 



2004 
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Schnittman fc Rezzolla II2005I ) 



Resonances in t 



Abramowicz et al. 



le ge odesic motion of a single particle have been investigated by 



(120031 ). The specific coupling force between radial and vertical 
oscillation was left unspecified. Their results for non-linear resonance were applied to 
accreting neutron stars. 

In this paper, we describe a coupling between spiral density waves in the accretion 
disk and epicyclic motions of test particles. It is divided in two parts. In Sec. [2] we specify 
the perturbation pattern used in our model. Then the equation of motions due to the 
perturbation are derived leading to some resonance conditions. In Sec. [3l the results are 
applied to GRS 1915+105 for which several components in the Fourier time analysis are 
predicted in agreement with the HF-QPOs detected for this object. We also put some 
constrain on the mass-spin relation for several BHCs. 



2. THE MODEL 

In this section, we describe the main features of the model, starting with a simple 
treatment of the accretion disk, assumed to be made of non interacting single particles 
orbiting in the equatorial plane of the black hole as was already done in the previous 
studies (see paper I and II). We again neglect the hydrodynamical aspects of the disk such 



as pressure. However, a detailed (magneto-)hydrodynamical treatment o 
of th e disk to gravitational or magnetic perturbations has been given in 



the r e sponse 



Petri 



fl2005c . 



20061 ). For a first approach to the problem, we neglect this refinement. Particles evolve 
in the unperturbed stationary background gravitational field of the black hole. In this 
simplistic approach, we want to point out the effect of a spiral density wave propagating 
radially outwards in the accretion disk. In order to account for general-relativistic effects. 
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we use the Kerr metric to perform our calculations. We will show that a local outgoing 
one-armed m-mode will induce a beat with the test particles and excite vertical oscillations 



parametrically. Indeed, 



MaoetaL 



(120081 ) showed that oscillatory modes are excited 
close to the innermost stable circular orbit (ISCO) and propagate outwards to several 
tens of Schwarzschild radii with constant frequency (nearly equal to the maximum radial 
epicyclic frequency) and without attenuation. They also found that these oscillations 
are locally super-Keplerian. Nevertheless, they studied only axisymmetric modes in a 
pseudo-Newtonian gravitational field. A more detailed study including full general relativity 
and asymmetric modes propagation would help to give precise quantitative results. In this 
paper, we will assume that such a s ymme tric modes exist and exhibit the same properties 



as the modes found in 



MaoetaL 



fl2008h . 



2.1. Equation of motion for a test particle 

We use the same procedure as the one described in paper I. However, in order to take 
into account the spiral structure of the density wave of azimuthal mode number m and 
frequency we have to change the equation of motion Eq. (7) in paper I by modifying 
the terms "cos(m {Q — f2„) t)" in such a way that the spatio-temporal dependency of 
the argument of "cos" implies a propagation of a pattern of azimuthal mode m at the 
sound speed Cg. We remember that Q is the local orbital frequency in the accretion 
disk. Therefore, in cylindrical coordinates {r,(p,z), we choose a phase dependency like 
ip = rmp — {t — r/cs), where t is the time. For a particle in circular orbit at a frequency 
fl, we have 

ip = mflt — (t — r /cs) (1) 

This means that viewed from the orbiting particle frame of reference, the perturbation is 
rotating at the speed Qp = \mQ — We find again a harmonic oscillator with a periodic 
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variation in the eigenfrequency of the system, driven by an external periodic force generated 
by the perturbation. The modulation being sinusoidal, the Hill equation specializes to 
the Mathieu equation, a well known ordinary d ifferential equation extensively studied in 



mathematical physics (IMorse fc Feshbach 



195 



2.2. Resonance conditions 

The rotation of the asymmetric pattern induces a sinusoidal variation of the vertical 
epicyclic frequency leading to the well known Mathieu equation for a given azimuthal 
mode m. To take into account general-relativistic effects, we use the characteristic orbital 
and epicyclic frequencies for the Kerr spacetime geometry, namely 



the orbital frequency Q 



r3/2+agjj GMbh r3/2+a 



the radial epicyclic frequency K^ = Vt\l\ — ^ + S — 3 ^ 



the vertical epicyclic frequency = -^/l — 4 + 3 || 



We introduced the mass of the black hole Mbh, its gravitational radius Rg = GMbh/c^, the 
adimensionalized radius and angular momentum, respectively f = r/Rg and a = a-Bn/Rg- 
We expect a parametric resonance related to this time-varying vertical epicyclic frequency. 
The resonance condition derived from the vertical equation of motion is 

\mn-n„\=2— (2) 
n 

where n > 1 is a natural integer. Note that the speed of propagation Cg in the term 
"(if: — r/cg)" drops out from the above resonance condition because it is a local description 
at fixed radius r. 
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In the non-rotating black hole case, a = 0, we have Kz = Therefore the parametric 
resonance conditions Eq. ([2]) splits into two cases, depending on the absolute sign 

n = — — (3) 

mn ± 2 

It differs by a factor 1/m from the resonance condition derived in paper I and II. This 
discrepancy comes from the choice of the phase dependency here being ip and describing 
a radially outward propagating pattern. There are some other resonance mechanisms 
happening in the accretion disk which are not of interest in the present study. Nevertheless 
we remind them for completeness. First, forced oscillations occur whenever the driven 
frequency is equal to the free oscillation frequency 

\mQ — = (4) 

They are a special case of the above mentioned parametric resonance, corresponding to 
n = 2. Second, for the corotation resonance to happen, we should have 

n = n^ (5) 

As a consequence, we find again the same resonance mechanisms as in paper I and II. 
Eq. ([3]) is also a good approximation for slowly rotating black holes (a ^ 1) because in this 

We now discuss in detail the case 5 = because it is analytically tractable and gives 
good insight on the resonance frequency even in the fast rotating case. Moreover, we 
checked that in case of a maximally rotating Kerr black hole, a = 1, the frequencies remains 
nearly the same as in the case a = because the resonances occur in regions far away from 
the ISCO, such that in these regions of the accretion disk, n^, ^ ^ whatever a. 

The predicted QPO frequency ratio Q/Q^ is shown in Table [1] for the first three 
azimuthal numbers m = 1, 2, 3 and the first four integers n. The highest orbital frequencies 
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where resonance occurs correspond to m = 1 and for the minus sign in Table [H They are 
equal to and 2Q^. Therefore the ratio of the strongest oscillations are in the ratio 3:2 
as observed in several black hole binaries. This particular ratio is not a direct consequence 
of the motion in general relativity but rather an intrinsic property of the parametric 
resonance. General relativity is only needed in order to describe this resonance mechanism 
correctly in the strong gravity regime. The third strongest resonance occurs when Q = 
overlapping the corotation resonance. However, if the low azimuthal numbers possess the 
strongest amplitude in the perturbation power spectrum, the m = 2 and m = 3 resonances 
are expected to be weaker. The most significant QPO frequencies are therefore ordered in 
the series 3:2:1, the "1" only appearing as a weak feature compared to the 3:2 HF-QPOs 
pair. A quantitative application of the model follows in the next section. 



3. DISCUSSION 



3.1. Application to GRS 1915+105 



We test our model with the well studied black hole GRS 1915+105 for which numerous 
observations of QPO features are available. Several types of QPOs can be identified in this 
binary system; low frequency QPOs (LF-QPOs) ranging roughly from 1 to 10 Hz, HF-QPOs 



for frequencies 



arger than roughly 70 



less than 1 Hz (IFender fc Belloni 



and z/2 = 165 Hz 



iz and very low frequency QPOs with frequencies 



Remillard et al. 



this object. For instance. 



2004f). HF-OPO frequencies are detected at z^i = 113 Hz 



Strohmayer 



2006). Ot 



ler LF-QPOs have also been detected for 
(j2001bl ) re ported a 40 Hz OPO seen sometimes 



simultaneously with the more stable 67 Hz QPO (IMorgan et al. 



observed with th e following increasing frequencies 27-40-56-67 Hz ( 



vanderKlis 



19971). OPOs have be en 



Belloni et al. 



2001 



intensively by 



2004f). Some LF-QPOs var ying between 1 and 15 Hz have been investigated 



Markwardt et al. 



(Il999l ). Their properties seem to be different from the 
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HF-QPOs, probably implying a different physical mechanism at work. We emphasize that 
it is not the scope of this paper to explain all the observed QPOs but just to explain or 
predict the HF-QPOs. LF-QPOs should be related to Lense-Thirring precession or to some 
other (magneto-)hydrodynamical modes in the accretion disk. 

Nevertheless, from the frequencies of the twin HF-QPOs, the constant angular pattern 
speed of the density wave is derived from our model by z/„ = VL^/2ti = v^jl = v^j'i ~ 56 Hz. 
Putting this value of the gravitational field pattern speed into Table [1] we get the results 
shown in Table O There is a problem on selecting the relevant and meaningful frequencies in 
this table. Indeed, the spiral density w ave travels outwards to several tens of Schwarzschild 



radii. According to 



MaoetaL 



(120081 ). the precise extension of the propagation of the wave 
with constant frequency and without attenuation depends on the accretion rate and on the 
viscosity (at least for their axisymmetric hydrodynamical modes). Assuming that the wave 
can travel only up to a radius Tout? for a black hole of mass Mbh and angular momentum 
aBH, the lowest excited frequency is 

_ Oio^ _ 32310 Hz 1 

2 TT Mbh/Mo aBH//?g + (rout/i?g)^/2 ^ ^ 

To give some orders of magnitude, let's use Mbh = lOM0,aBH = O,rout = 30i?g, we find 
^'low = 19.7 Hz. Which frequencies to select depends strongly on the spiral wave properties, 
how far they can propagate, which azimuthal numbers m are excited, at which amplitude 
and so on. The detailed investigation is left for future work. We can only claim that the 
observed QPOs are retrieved by our model. How to select them remains quantitatively 
unclear. 

Let's discuss however the predicted QPO frequencies. The first two strongest resonances 
for each azimuthal mode m with the plus sign used in Eq. ([3]) match with good accuracy 
the frequencies observed in GRS 1915+105. We find namely, in decreasing order, 168.0, 
112.0, 56.0, 42.0 and 28.0 Hz for the most significant values. Some other resonances should 
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also appear, but due to their higher order (higher integer n) and location farther away from 
the inner edge of the disk (thus an amplitude of the excitation decreasing with radius), they 
do not possess a significant growth rate relevant for the study presented here. Moreover, 
X-ray emission from the outer part of the disk is fainter, thus more difficult to detect. As 
a consequence, the parametric resonance induced by a spiral density wave passing through 
an accretion disk predicts the HF-QPOs as well as some LF-QPOs (as a byproduct) in a 
single unified picture. However, there is no way to predict the angular momentum of the 
black hole without some knowledge of its mass. We can derived a relation mass-spin to 
constrain th e BH parameters, see next section. Nevertheless, it is worthwhile noting that 



Kato 



(12004 ) was lead to the conclusion that this microquasar is well described by a non 
rotating geometry indicating that the black hole angular momentum is weak, obh ^ Rg by 
fitting the HF-QPOs and several LF-QPOs with a resonant interaction between non linear 
oscillations and warp modes in the accretion disk. He found a^n = — 0.15 i?g. 

Note also that the LF-QPOs (27-41-56 Hz) are more difficult to detect than the 
HF-QPOs (113-165 Hz) because they are located in regions of the accretion disk where 
emission is fainter. They correspond also to higher azimuthal modes (m = 2, 3 compared 
to m = 1). Assuming that the main wave possesses an m = 1 structure, the amplitude of 
the LF-QPOs is expected to be small compared to those of the HF-QPOs. Furthermore, 
it seems that the presence of these features depend on the emission state of the black 
hole. We would then expect the accretion rate to have an influence on the efficiency of 
the parametric resonance occurring in the system. However, the constancy of the QPOs 
frequency strongly supports the fact the they are closely related to the black hole properties 
and independent of the flow in the surrounding accretion disk. A hydrodynamical general 
relativistic description is required to verify this assessment and to check the behaviour of 
the accretion r ate on the p r esence of certain QPOs. This would be the continuation of the 



work begun by 



MaoetaL 



(120081) 
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3.2. Mass-spin relation for BHCs 

Strictly speaking, at this stage, our model cannot estimate independently the angular 
momentum and the mass of the black hole. Nevertheless, applying it to some BHCs, we are 
able to give some constrain on their mass, assuming a non-rotating or a maximally rotating 
black hole. 



Following the work done by iMao et al. I (120081 ). we assume that the spiral density 



wave is launched close to the ISCO, at the location where the radial epicyclic frequency 
is maximal, note this radius rmax- The frequency of the spiral density wave is equal to 
the local orbital frequency at rmax, thus the speed of the density perturbation pattern is 
^('^max, c^bh)- Kuowiug the fundamental frequency fl„ in the accretion disk for a given 
black hole, the resonance condition Eq. ([2]) puts constrains on the relation between mass 
and angular momentum by imposing fi(rmax, o-bh) = 

This mass-angular momentum relation is shown in Fig.[T]for four BHCs (GRS1915+105, 
GRO J1655-40, XTE J1550-564, H1743-322) for which the fundamental frequency is known. 
If we assume a black hole weakly or mildly rotating, i.e. obh ^ 0.5 G M-Qn/c^, we find that 
for GRS 1915+105 its mass satisfies 13 Mq < Mbh < 20 Mq. The same model applied to 
other BHCs gives for GRO J1655-40 a mass 5Mq < Mbh < 7 Mq and for XTE J1550-564 
a mass 8 Mq < Mbh < H Mq. This is consistent with other independent estimations of the 
BH mass. Finally for H1743-322, we found the following bounds, 9 Mq < Mbh < 13M0. 



4. CONCLUSION 

In this paper, the consequences of an outgoing spiral density wave on the evolution of 
a test particle orbiting around a Kerr black hole have been explored. The twin peaks ratio 
around 3:2 for the kHz-QPOs is naturally explained by the parametric resonance model. 
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0.2 0.4 0.6 0.8 1 

Fig. 1. — The relation between mass Mbh (normalized to the mass of the Sun Mq) and 
angular momentum aeH for different fundamental excitation frequencies corresponding to the 
four BHCs, namely, GRS1915+105 with 56 Hz, GRO J1655-40 with 150 Hz, XTE J1550-564 
with 92 Hz and H1743-322 with 82 Hz. 
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The connection to lower frequency QPOs is also clearly demonstrated. It is an extension 
to accreting black hole for the model already suggested in neutron star and white dwarf 
binaries. 

From the analysis of HF-QPO s, we are also able t o constrain the mass and the spin 



of the hole as already suggested by 



Torok et al. 



(12005h . In the case of GRS 1915+105, it 



appears that the 56 Hz feature in the Fourier time analysis is a fundamental frequency of 
the black hole from which the other QPOs can be derived. This feature should therefore be 
related to the intrinsic properties of the black, namely, its mass and its angular momentum. 
The way to associate with and Mbh remains unclear and needs further investigation 
but spiral density wave excitation cl ose to the I S CO a t nearly the maximum of the radial 



epicyclic frequency as suggested by 
further investigations. 



Mao et al. 



(I2OO8I ) is a interesting idea that needs 
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Table 1: Ratio of orbital frequency to the perturbation angular pattern speed at the radii 
where the parametric resonance is expected to occur. Left part of the table corresponds to 
the minus sign whereas the right part corresponds to the plus sign taken in Eq. ([3]). 



m 


n 


(-) 






n 


(+) 






1 2 


3 


4 


1 


2 


3 


4 


1 


-1 — 


3 


2 


1/3 


1/2 


3/5 


2/3 


2 


— 1 


3/4 


2/3 


1/4 


1/3 


3/8 


2/5 


3 


1 1/2 


3/7 


2/5 


1/5 


1/4 


3/11 


2/7 
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Table 2: QPO frequencies predicted for GRS 1915+105 from the orbital motion wiiere para- 
metric resonance occurs. The most interesting frequencies are marked in bold face. 



i/QPO (in Hz) for GRS 1915+105 



m 




n 


(-) 




n(+) 




1 


2 


3 


4 


1 


2 


3 


4 


1 


-56.0 




168.0 


112.0 


18.7 


28.0 


33.6 


37.3 


2 




56.0 


42.0 


37.3 


14.0 


18.7 


21.0 


22.4 


3 


56.0 


28.0 


24.0 


22.4 


11.2 


14.0 


15.3 


16.0 



